Flavour symmetry of mesonic decay couplings 



Eef van Beveren 

Departamento de Fisica, Universidade de Coimbra 
P-3000 Coimbra, Portugal 
eef@malaposta.fis.uc.pt 

George Rupp 

Centra de Fisica das Interacgdes Fundamentais 
Instituto Superior Tecnico, Edificio Ciencia 
P-1096 Lisboa Codex, Portugal 
george@ajax.ist.utl.pt 

PACS number(s): 14.40. Cs, 12.39.Pn, 13.75.Lb 



hep-ph/980624g 



February 1, 2008 



Abstract 

We present flavour-symmetric results for the couplings of quark-antiquark systems to 
meson-meson channels in the harmonic-oscillator expansion. We tabulate their values for 
all possible open and closed decay channels of pseudo-scalar, vector and scalar mesons. We 
compare the predictions of a model that employs these flavour-symmetric couplings, both 
with the results of a model which uses explicitly flavour-dependent couplings, and with 
experiment. 



1 Introduction 



Particle interactions are described by point-particle vertices in fundamental theories. Quarks, the 
basic particles for strong interactions, are point objects, to our best knowledge, hence assumed 
to interact via point-particle vertices in the existing theories: through a quark-gluon vertex in 
Quantum Chromodynamics (QCD) through a four-quark vertex in the Nambu-Jona-Lasinio 
model (NJL) §. 

QCD exhibits good agreement with experiment, qualitatively for low and medium energies, and 
moreover quantitatively at high energies , whereas NJL shows good agreement with experiment 
only for energies below 1 GeV [Q. So in the energy interval crucial to meson physics, i.e., ranging 
from the two-pion threshold to energies as high as the states in the bottomonium system, no 
fundamental theory possesses a satisfactory descriptive power: not QCD, because the relevant 
momentum transfers are too low and thus the effective color coupling constant is too large 
for a perturbative approach, not NJL, because the energies are too high. Consequently, for a 
quantitative description of the spectra and scattering of mesons and baryons, neither of the two 
theories has sufficient predictive power for the time being. Therefore, the use of quark models is 
still opportune in this domain of hadronic physics. 

Now, ideally a quark model should be derived from QCD, but this is rather Utopian as yet. As a 
matter of fact, not even a direct relation between QCD and confinement has been established and 
so confinement usually has to be imposed on the valence quarks of the model ^ . Different models 
follow distinct strategies to achieve this and the manner in which confinement is approached 
distinguishes models among each other. Moreover, each model has its own very specific purpose, 
often not mentioned in too much detail by the authors, which makes it difficult to compare 
models. For instance, there exist heavy-quark potential models made to measure in order to 
reproduce, with great accuracy, the radial and angular spectra of charmonium and bottomonium, 
as well as the electromagnetic properties of these systems. But if the same potentials are used 
in the light quark sector, the results are normally quite bad, especially for radial excitations, 
and even possible relativistic corrections are insufficient to cure the discrepancies. Conversely, 
sophisticated relativistic models for the light mesons usually fail to reproduce the correct radial 
spacings in the charmonium and bottomonium spectra. 

On the other hand, most quark models treat hadrons as manifestly stable bound states of 
quarks, simply ignoring the fact that most hadrons are resonances, some of them even extremely 
broad, so as to make their very existence questionable. The standard justification is the con- 
jecture that the effect of strong decay will be to produce predominantly imaginary mass shifts, 
thus allowing to first fit the real parts of the spectra and then to treat the hadronic widths a 



posteriori, with perturbative methods. However, we know from fundamental principles in scat- 
tering theory that real shifts are generally of the same order as or even larger than the imaginary 
ones. Moreover, hadronic loops, i.e., virtual decay channels, give rise to attractive forces, so that 
the shifts due to, in principle, all closed decay channels, must be added up so as to produce a 
negative mass shift. So not even the true bound states can be treated as pure quark states. The 
usual excuse is the unsupported assertion that the effect of closed channels will be negligible, 
except near threshold. 

However, the Nijmegen unitarised meson model (NUMM), devised to simultaneously describe 
meson spectra and meson- meson scattering, from the light pseudoscalars and vectors , via the 
usually awkward scalars 0, all the way up to the bb sector showed that both premises are 
indeed wrong: real shifts are generally comparable with or larger than the imaginary ones, and 
the damping of closed channels is insufficient to make their influence on the ground states of the 
spectra negligible. On the contrary, due to the nodal structure of the radial wave functions and 
the mentioned additivity property, the shifts — real and negative - of ground states are usually 
largest . Furthermore, no drastic enhancement takes place near threshold, so that these states 
cannot be singled out [|T0[] . 

Having come to the conclusion that, for a truly quantitative description of mesonic spectra, 
one must include the coupling to meson-meson channels, the crucial questions to be raised are 
how to calculate the involved coupling constants and which two-meson channels to take along. 
Here, one should step back and have another look at the QCD Lagrangian. Realising that, at 
least qualitatively, there should be no obvious disagreement between QCD and whatever meson 
model to be used, we are led to respect manifest flavour blindness. This will impose stringent 
conditions on how couplings can be computed and how to select classes of decay channels, since 
obviously one cannot take into account an infinite number. 

Models which describe the scattering of mesons (and/or baryons) often imitate the funda- 
mental theories in the sense that interactions take place via effective point-particle vertices. 
However, meson (and baryons) are composite systems, built out of strongly interacting valence 
quarks, glue, and a quark-antiquark sea. So it seems obvious that, when mesons (and baryons) 
are considered point objects, some information must get lost. In this paper we will demonstrate 
that this can indeed be the case and how it manifests itself in the flavour non-independence of 
the such-described strong interactions, when not dealt with carefully. 

Point interactions are a powerful tool in constructing theories that not only consider relativistic 
kinematics, but also take into account the property of particle creation and annihilation. However, 
in applying point interactions to composite systems, one should include all hidden degrees of 



freedom. Flavour is just one such degree of freedom. Angular momentum and spin are others 
which should be properly included. At present, it is opportune to model the internal degrees 
of freedom and next to integrate them out for the determination of effective point couplings. A 
consistent way of doing so, which moreover preserves flavour independence, is described below. 

The organisation of this paper is as follows. In section ^, we discuss the general philosophy 
behind a simple model for flavour symmetry. This model is then exposed in section |^. The 
intensities of the three-meson vertices for meson decay into meson-meson pairs are given in section 
^. Results are discussed in section ^ The consequences of flavour (in) dependence are studied for 
two different, though similar, models in section ^ Some essential formulae are collected in the 
appendices and |B[ 

2 Flavour symmetry 

Since strong interactions are independent of flavour, the probability to create a quark- ant iquark 
pair out of the vacuum cannot depend on the flavour of the quark and the antiquark. However, 
it obviously depends on the masses involved. But if for a moment we assume that the flavour 
masses, or at least the effective quark masses in the relevant energy interval of the three lowest 
flavours, up, down and strange, are equal, then the corresponding probabilities of pair creation 
should be equal. Let us apply this principle to the strong coupling of a meson to a pair of 
mesons. Here, we assume that the related strong decay processes are triggered by the creation 
of a flavourless quark-antiquark pair. 

In order to set the picture, we consider a simple model in which the initial meson is described 
by a confined quark-antiquark system of any flavour, given by 

ab , (1) 

where a and b represent any of the three flavours under discussion, and the final pair of decay 
products is described by a system of two freely moving mesons, which represent any of the three 
combinations 

(au) + iub) , (ad) + (db) , (as) + (sb) . (2) 

Also, let us, for a moment, assume that no further quantum numbers are involved. Then, 
flavour symmetry demands that the probabilities are equal for the ab system to decay into any 
of the three channels of formula (|^). In particular, when under full flavour symmetry the six 



mesons represented in formula have all the same mass, then the experimental results for the 
decay of system (|I|) into any of the three channels of formula (Q) should be indistinguishable. 

Let the three decay coupling constants of the process under consideration be represented by, 
respectively, 

g{a,b;u) , g{a,b;d) , g{a,b;s) . (3) 
Then, assuming full flavour symmetry, we have the identities 

{g{a,b;u)}'^ = {g{a,b; d)}"^ = {g{a,b; s)}"^ , (4) 
and, moreover, for the total decay intensity r{a,b) the relation 

T{a,b) = A [{g{a,b;u)}'^ + {g{a,b;d)Y + {g{a,b; s)}'^\ , (5) 

where the proportionality factor A is also completely flavour symmetric, which means constant 
here. Furthermore, one has that, under flavour symmetry, T{a, b) must be independent of the 
flavours a and b. 

Unfortunately, quarks are fermions and mesons are spatially extended systems, and hence 
spin and spatial quantum numbers do play an important role in the decay of a meson into a 
pair of mesons. Nevertheless, it remains possible to construct coupling constants which have the 
property that the total decay probability is independent of the flavour of the decaying meson in 
the limit of equal masses, as we will see below. 



3 Modelling full flavour symmetry 



When a normalised wave function ifj is expanded on a complete orthonormal basis for n = 
0, 1,2,.. ., according to 

oo 

i' = ^^^^ ' (6) 

n=0 

one has for the expansion coefficients c„ the property 

oo 

E ICnl' = 1 ■ (7) 
n=0 

It is exactly property (|^) that leads to flavour symmetry. 

Let us consider a system of two quarks and two antiquarks, like any of the three combinations 
of formula (^. One complete basis for the Hilbert space of such a system can be constructed 
by taking products of the internal wave function of the ab system, specifying thereby its internal 
spatial and flavour quantum numbers, the internal wave function of the qq (either uu, dd or ss) 
system, and the relative wave function of the two subsystems. Another complete basis for this 
Hilbert space consists of products of the internal wave function of the aq system, the internal 
wave function of the qb system, and the relative wave function of those two subsystems. Any 
wave function describing one of the two-quark-two-antiquark systems (0) can be expanded in 
either of the above-defined bases. 

Such an expansion takes a particularly manageable form if the four partons are supposed 
to move in a harmonic oscillator potential with universal frequency. In that case the spatial 
quantum numbers are linearly related to the total energy of the system, which gives rise to finite 
bases at each energy level and hence to finite expansions. The fiavour-symmetry condition (|^) 
then becomes a finite sum, which makes it easy for verification. Furthermore, the restriction to 
harmonic oscillators is not a real limitation, since any other basis can always be expanded in the 
corresponding harmonic oscillator basis, {n}, according to 

{MiM2\V\M) = ^ (M1M2 (n'|K|n) (n|M) . (8) 

Here, V represents the interaction Hamiltonian which describes the transitions between the quark- 
antiquark system ab and the two-meson channels. We assume that the spatial, or momentum- 
dependent, part of the matrix elements {n'\V\n) is fiavour independent and that the fiavour- 
dependent parts are constants. 

The expansion of a particular many-particle wave function into a specific basis for well-defined 
subsystems, or recoupling, has been studied a great deal in the past. The related coefficients for 



the harmonic oscillator basis are known as Moshinsky brackets. Moshinsky brackets are well- 



known coefficients of recoupling in Nuclear Physics; see Ref. |[TT]] for their definition. The group- 
theoretical implications of parton recoupling in the harmonic-oscillator approximation have been 



studied exhaustively in Ref. [|12[. Their application to meson decay has for the ffist time been 
formulated in Ref. |jl3|. A full generalisation for the spatial part of the recoupling constants, 
which includes all possible quantum numbers for any number of (bosonic) partons, can be found 



in Ref. [H]. The inclusion of fermionic and flavour degrees of freedom, which leads to an analytic 
expression for the coupling constants of any meson to any of its two-meson real or virtual decay 
channels, is given in Ref. |T5[ , for the case that the new valence pair is created with ^Pq quantum 
numbers. A Fortran source program is available on request. 



4 Coupling constants for three-meson vertices 



Within the above-outhned formahsm, we assume that mesons can be classified by the quantum 
numbers of their valence constituent quark- ant iquark distributions, i.e., 

meson {j,M,i,s,n,M) . (9) 

The quantum numbers j and M in formula (H) represent, respectively, the spin of the meson and 
its z-component. Alternatively, j represents the total angular momentum of the relative motion 
in the quark+antiquark system which describes the meson. The quantum numbers £, s and n 
stand, respectively, for the orbital angular momentum, the spin, and the radial excitation of the 
constituents of the meson. Finally, Ad represents the 3x3 flavour matrix which indicates the 
valencies of the quark and the antiquark. 

Here, we study the decay intensities for the following processes: 



meson{J,J^,i,s,n,Mc) — ^ meson (ji. Mi, £i, Si, rii, A^a) + meson (j2, M2, £2, ^2, n2, A^s) . 

(10) 

This is not a completely satisfactory notation, since the spin z-components Mi and M2 of the 
decay products in formula (|TUp are not supposed to be observable and, moreover, the quantum 
numbers which characterise the relative motion of the decay products are not specifled in formula 
(|10D, despite being equally important. Let us indicate the orbital quantum numbers of the two- 
meson system by means of an index, r, and hence denote the orbital angular momentum of 
the two mesons by ir, the total spin by Sr, and the radial excitation of the relative motion in 
the two-meson system by Ur. The total angular momentum of the two-meson system and its 
^-component are, due to angular- momentum conservation, given by J and Jz, respectively. 

The decay probability for the process (|10D is then, following the formalism developed in 
Ref. [jl5|, given by the following matrix element: 



( J, Jz, jlJl, Si, ni,j2, 4, S2, 4, Sr, Ur, A,B\J, J^, i,s,n,C) = (11) 
= Tr{ MaMbMc^ {J,Jz,jl,il,Si,ni,j2,i2,S2,n2,ir,Sr,nr\J,Jz,i,S,n,a.ABc) + 

+ MbMaMc^ {J,Jz,ji,ii,Si,ni,j2,i2,S2,n2,ir,Sr,nr \ J, J z, s,n, olbac) } • 
The spatial parts in each of the two terms of the transition element (pTf) are denoted by 



{J,Jz,jl,il,Si,ni,j2j2,S2,n2,ir,Sr,nr\J,Jz,(^,S,n,a) , (12) 

and are defined and explained in Appendix |^. 

Since for many purposes it is sufficient to liave flavour-symmetric coupling constants for the 
corresponding strong decay channels of pseudoscalar, vector, and scalar mesons, we tabulate the 
probabilities for the three-meson vertices of those decay processes; for a pseudoscalar meson in 
Table ^ for a vector meson in Table |^, and for a scalar meson in Table ^. In order to maintain the 
tables as condensed as possible, we represent mesons by symbols and by their quantum numbers. 
Since we assume that isospin is indeed a perfect symmetry, we may represent all members of an 
isomultiplet by the same symbol, for which we just have chosen the letters and numbers t, d, 8 
and 1, according to the identification given in Table |1]. Now let us just analyse one horizontal line 
of one of the three tables, to make sure that the reader understands what the numbers represent. 
Let us take the fourth line of Table |^. In the first column we find four zeroes, representing the 
internal spatial quantum numbers j, i, s, and n of the first decay product. Mi, which hence 
characterises a meson out of the lowest-lying {n = 0) pseudoscalar nonet. In the second column 
we find similarly that the second decay product, M2, represents a meson out of the lowest-lying 
vector nonet. In the third column we find the quantum numbers for the relative motion of 
Ml and M2, i.e., P-wave {ir = 1) with total spin one (s^ = 1) in the lowest radial excitation 
{rir = 0). Since the table refers to the real or virtual decays of the lowest-lying pseudoscalar 
meson nonet {Jisn = 0000, indicated in the top of the table), the next four columns refer to its 
isotriplet member, which is the pion. We then find that the pion couples with a strength 
to the tt (isotriplet-isotriplet) channel, which, following Table |I] and the above-discussed particle 
assignments to Mi and M2, i.e., pseudoscalar and vector respectively, represents in this case the 
Tip channel. Following a similar reasoning, we find that the pion couples with a strength 
to KK* . The total coupling of a pion to pseudoscalar- vector channels is given in the column 
under T by -y/l/4, which is the square root of the quadratic sum of the two previous couplings, 
i.e., ^1/6+1/12. 

The next set of coupling constants refer to the real or virtual (actually only virtual) decays of 
a kaon. We find y'l/S to trf, which represents both of the possibilities pseudoscalar (isotriplet) 
+ vector (isodoublet), i.e., tiK*, and pseudoscalar (isodoublet) + vector (isotriplet), i.e., Kp, 
each with half of the intensity that is given in the table, and therefore one has for the kaon the 
coupling constants to nK* and to Kp. Next, we find in the table that the kaon 

couples with to (iS, which represents both of the possibilities, pseudoscalar (isodoublet) + 

vector (Sf/s-octet isoscalar), i.e., K + some admixture of the uj and mesons, and pseudoscalar 



(S'?73-octet isoscalar) + vector (isodoublet), i.e., some admixture of ?7's + K* , each with half of 
the intensity that is given in the table, and hence one extracts for the kaon the coupling constants 
y^l/16 to + (tu, (p) and ^1/16 to (r/, rj') + K* . The kaon does not couple to the dl channels 
in pseudoscalar + vector, which represent the channels with one isodoublet and one S'?73-singlet. 
The total coupling for the kaon to its pseudoscalar + vector decay channels sums up to ^1/4, 
as one verifies in the column under T. The next two sets of coupling constants refer similarly to 
the decay modes of the isoscalar, either S'L^s-octet or S'f/3-singlet, partners of the pseudoscalar 
nonet. Mixings can be done by hand as examplified in Appendix |B[ 

Notice that for the SU ^-ociei members one has flavour symmetry for each horizontal line in 
the tables. This does not go through for the S'f/3-singlet partners, with the exception of scalar 
meson decay (Table where all horizontal lines have the same total coupling in each subsection 
of the table. However, all columns under T sum up to 1, representing full flavour symmetry once 
all possible decay channels have been accounted for. 



5 Results 

The pole structures of the scattering matrices for P-wave meson-meson scattering, or equiva- 
lently, the radial spectra of heavy and light pseudoscalar and vector mesons, were studied in the, 
largely non-relativistic, coupled-channel NUMM, published in Ref. hereafter referred to as 
B83, in which the authors parametrised confinement by a universal frequency, the same for all 
flavours, including charm and bottom. The universal frequency and a flavour-independent overall 
coupling constant, representing the probability for the creation of a ^Pq light quark- ant iquark 
pair, were sufficient to obtain theoretical predictions for phase shifts and scattering cross sections, 
or equivalently, for central resonance positions and widths, which were in reasonable agreement 
with the data. All relative couplings were exactly taken as given in Tables ^ and ^, though some 
of these had been derived in a more empiric way, and then extended in order to also include 
the heavy-quark systems. This extension is quite trivial and will not be discussed here. The 
only flavour non-invariance came from the quark masses and the two-meson thresholds, all other 
ingredients were the same for all flavours. Of course, many of the decay channels were omitted, 
assuming their thresholds to be high enough in energy, so as not to have too much importance 
for the details of the scattering processes at much lower energies. But this is only a practical 
ingredient, not to be confused with flavour breaking. 



In Ref. |]T6[, the electromagnetic transitions in the charmonium and bottomonium systems 
were studied, using the quark and meson distributions from B83, with good results, indicating 
that not only the pole structures of the scattering matrices, but also the related wave functions 



stood the confrontation with experiment. 

In Ref. [0], hereafter referred to as B86, the pole structure of the scattering matrix was 
inspected for 5- wave meson-meson scattering. Since the model was the same as for Ref. 
using exactly the same universal frequency, flavour- independent overall coupling constant, and 
quark masses, the calculated phase shifts and scattering cross sections could be considered genuine 
theoretical predictions. The agreement with the data was unexpected, especially because it had 
not been the objective of the model, neither was the model constructed towards fitting the S'-wave 
scattering data. All relative couplings were exactly taken as given in Table ^. Also here, only 
those two-meson channels were taken into account which contain members of the lowest-lying 
pseudoscalar and vector nonets. That such a procedure does not break flavour invariance may 
be explicitly verified by checking the first, third, and last line in Table ^. 

The observed flavour independence of strong interactions is a very important ingredient for 
low-energy hadron physics and woven into the NUMM; first, by the universal frequency, which 
makes the ratio of the kinetic term and the potential term of the model flavour independent 
and hence also the level splittings; second, by the intensities of the three-meson vertices for the 
coupling to the various decay channels. 



6 Comparison of two models 



As stated in the introduction, it is not easy to compare meson models, but here we will pay 



attention to the comparison of the NUMM with a model hereafter referred to as T95, which 
is tailor-made for scalar mesons or, in other words, for S-wave meson-meson scattering. The 
latter model, a revised version of the Helsinki unitarised quark model (HUQM), was confronted 



with experiment in an analysis published in Ref. [T^, hereafter referred to as TR96. Based on 



the good agreement of its theoretical predictions with the available experimental phase shifts, 
one may be inclined to accept all further conclusions presented in the same publication, such as 
the existence and location of resonances. However, the authors failed to find the complex-energy 
pole corresponding to the established /o(1500) resonance. Furthermore, they also did not find a 
light Kq, i.e., the old k, somewhere between 700 and 1100 MeV. Although the latter resonance 
is not (yet) established experimentally, it has recently received renewed phenomenological and 
theoretical support |1^, ^ ^ (see also Ref. ||22[). Moreover, its absence in Nature would imply 



a breaking of the conventional nonet pattern for mesons. On the other hand, if a light Kq is 
confirmed, then there exist unmistakable experimental candidates for two complete scalar nonets, 
as predicted by B86. So it is intriguing to figure out why model T95/TR96, which is very similar 
in its philosophy and also in observing a resonance doubling, at least for some states, does not 
reproduce this resonance. 

In Table ^, we collect the intensities for strong scalar-meson decay into a pair of pseudoscalar 
mesons, under the assumption that pions, kaons, and eta's have equal masses, as given by models 
B86 and T95/TR96. For the purpose of comparison, we have multiplied the values given in B86 
by a constant factor. The resulting values can also be read from the first line of Table ^ when 
renormalised (i.e., multiplied by a factor 24), and when isoscalar mixing has been dealt with as 
outlined in Appendix 0. Now notice that, by using formula (^, the total decay intensities stem- 
ming from Ref. [0 become equal to A for all scalar mesons, as demanded by flavour symmetry. 



r(ao) = r(/t) = r(/o,nn) = r{fo,ss) = A . (13) 

For model T95/TR96, the comparable intensities are derived from a point-particle approach 
to the three-meson vertex, which results in coupling constants given by 

Ti (MaMbMc) , (14) 
where A, B and C stand for the three mesons involved at the vertex C —>■ AB, and Aix is the 



3x3 flavour matrix for meson X. Flavour symmetry as from Eq. (JT^, and using formula 
yields in this case the flavour-dependent result 



r(ao) = r{K) = ^Tifonn) = {fo ss) = A , (15) 
5 4 



in contrast with the results shown in Eq. 

It might be surprising that such a seemingly flavour-symmetric vertex leads to flavour non- 
independence when applied to strong meson decay. However, if one does not take into account 
the internal structure of the various quark- ant iquark systems involved, then normalisations, so 
essential to wave functions or distributions, are swept under the rug, and thus the above vertex 



intensity (14) leads to flavour-dependent results. It appears to be for this reason that the authors 
of TR96, which base their calculations on the coupling constants from T95, do not observe any 
resonance doubling for the isodoublet and one of the two isoscalars, and therefore miss the 
-f^Q (700-1100) and /o(1500) poles needed to complete two scalar meson nonets. 

The normalisation factors that are relevant to formula ([T3|) are given in formulae ([l^) and ([T7| ) 
of Appendix 0. In Appendix we show how they lead exactly to the factors | and | which are 
necessary to compensate the flavour-dependence of formula (p!5|). 



A Rearrangement coefficients 



The spatial parts of the matrix elements ( ]TTD are, following the formalism developed in Ref . , 
in the approximation of equal flavour masses just given by Clebsch-Gordonary and some overlap 
integrals, amounting to 



{J, JzJl, 4, Si, ni, js, 4, S2, n2, ir, Sr, Ur | J, J^, i, s, n, olabc ) 
(J, Jz,jl, il, Si, ni, j2, ^2, S2, n2, 4, Sr-, Ur | J, J^, S, H, CXbAC ) 



E 
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(16) 



where the sum is over all /x's, m's, and M's that appear in the formula, where 



{C \SU (3)flavour-singlet ) 



for 


\C) = \uu) 


for 


\c) = 


for 





for \C) orthogonal to the S'f/(3)flavour-singlet state , 

ddj , or \ss) , 

I ^\uu) + ddj + |ss)| 



(17) 



and where 



;i8) 



The central part of formula (|I6]) is constituted by the rearrangement coefficients, which can 
be given by the following diagramatic representation 



/ n i m£ 
1m 
^000 



rii £i nil \ 



cx 




(19) 



-"33- 



The upper-left external line carries the relevant quantum numbers of the initial meson in formula 
(|IOD, the middle-left external line the quantum numbers of the ^Pq gg-pair, and the lower- left 
external line the quantum numbers of the relative motion of the two quark- ant iquark systems, 
which, to lowest order, is supposed to be in its ground state. The external lines on the right-hand 
side of the diagram carry the quantum numbers of the decay products of formula ( p!0[ ) and their 
relative motion. 

As explained in Ref . [l^ , each of the internal lines ij of diagram ([T9|) carries the set of quantum 
numbers {rijj, mjj}, over all possibilities of which must be summed, thereby respecting partial 
quantum- number conservation at each vertex, i.e.. 
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where the angular- moment a recouphng coefficients are defined by 
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and where the a-matrices, for the case of equal constituent fiavour masses, are given by 
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The allowed values for the quantum numbers of the internal lines of diagram (|T9|), given by riij 
and lij in formula (^), are non-negative integers and hence, because of partial quantum number 
conservation at each vertex of the diagram, which is moreover expressed by the Kronecker delta's 
in formula (|20|), we find 



ni2 



^22 = ^32 



^13 



^23 = ^^33 



'13 



'23 



^^33 







(23) 



which, also substituting the CK-matrices (]2^), simplifies the expression for the rearrangement 
coefficients to 
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Notice that the Kronecker deltas in formula i^M amount to 



J^Sl ^Z2 Ir 





(24) 



2n + i + l = 2{ni+n2 + nr) +ii+£2 + ir , (25) 
which is precisely the important relation that limits the number of possible decay channels. 

Moreover, for an initial pseudoscalar or vector meson out of the lowest-lying flavour nonets, 
one has in formula (|ll]) for the qq quantum numbers n and £ that 



n = i = . (26) 

Consequently, through the use of the Kronecker deltas in Eq. (^4]) , we find for the quantum 
numbers of the internal lines of diagram (p^) that 

nu = n2i = nsi = in = £21 = £31 = , (27) 

and moreover 



m = n2 = = and £1 + £2 + £r = I ■ (28) 

Relations (^) can be checked against the first three colums of Tables (H) and (^, where for 
all possible channels the radial excitations ni, 77.2, or Ur vanish, and, moreover, the sums of £1, 



^2, and ^r equal 1. This is a consequence of formula (|25| ) and drastically limits the number of 
possible quantum numbers and hence decay channels. 

For the relevant rearrangement coefficients we find, using formula (^), in this case 
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For the decay of a meson out of the lowest-lying scalar nonet, one has 

n = and i = s = 1 , 



(30) 



and hence, by the use of formula (p^ , one obtains for the relevant rearrangement coefficients in 
this case 
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B Mixing for scalar mesons 



Since we assume ^Pq quantum numbers for the creation of a qq pair out of the vacuum, the mixing 
of the isoscalar flavour-nonet members for lowest-lying scalar meson decay is not completely 
trivial. So we will outline here some of the necessary ingredients. 

In order to simplify the discussion, let us separate the normalisation factor and the summations 



in formula (0). Moreover, the two Kronecker deltas under the square root in formula (|T^) do 
not vanish for the lowest-lying scalar mesons. Therefore, let us denote 



(0, 0, ji, £i, si, ni, ja, 4, S2, ns, 4, Sr, rir |0, 0, 1, 1, 0, cxabc ) 

(0, 0, ji, ii, Si, ni, j2, 4, S2, n2, 4, Sr, rir |0, 0, 1, 1, 0, (XbAC ) 



(32) 
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1 + {C |S'f/(3)flavour-singlet) 




where (ABC) stands for the upper summation in Eq. ([TBD and (BAC) for the lower. 

As one notices from formulae (p!6D, (|19|), (^) or just only from formula (|3TD, the transition 
coefficients {ABC) and {BAC) do not , for full SU{3) flavour symmetry, i.e., for equal up, 
down and strange quark masses, depend on the flavour contents of the three mesons A, B, and 
C involved in the transition process (p^Of ), but just on the orbital and intrinsic spin quantum 



numbers of the system, which circumstance is also expressed by the notation of formula ([l^) . In 
fact, for the case of equal quark masses, those transition coefficients are equal for each different 
set of spatial quantum numbers, up to a sign. This sign is positive for all possible couplings in 
the case of the lowest-lying scalar mesons. Consequently, since for Table ^ only the transitions 
between scalar mesons and pairs of pseudoscalar mesons are relevant, and for mixing in general 
only the restriction to a specific set of spatial quantum numbers has to be considered, we may 
put here 



(ABC) = {BAC) . (33) 

For the transitions of the flavour-octet members to pairs of mesons, one has, according to 
formula ([T7|), a unity normalisation factor. Let us study then the matrix elements for a repre- 
sentant, ud, of the isotriplets, which is denoted by t in Table ^, coupled to an isoscalar, 0, and 
an isotriplet, for which we also take as a representant the ud state, and which, moreover, is also 
denoted by t in Table ^, i.e.. 



(udj 



ud 



If (f) represents the flavour-octet-member isoscalar ^g, we find 



ud 



ud) w I (uu + dd — 2ss 



ud 



Obviously, the matrix element for the ss contribution to 08 vanishes, hence 



(ud^ (ps ud"^ = y|^{((wc?^ (uu) ud^ + (^(ud'j {dd}j ud^^ 



which, also using Eq. (PBD, gives 
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(35) 



Following a similar reasoning when in Eq. (34) represents the flavour- singlet isoscalar 0i, 
we find for its matrix elements the result 



ud 



udj 



ud) J I (uu + dd + ss 



ud 
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(36) 



In Table ^, for the quadratic matrix elements under t8 and tl in the sector isotriplets, one may 
verify the factor 2 that follows from formulae ( P3| ) and (PB]), 

For the ideally- mixed isoscalars, (pn (non-strange) and 0^ (strange), defined by 
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one finds the matrix elements 



(ud) 0„ 
ud) 0s 



ud) = V2{ABC) = V3((mc/) 
ud) = . 



ud 



^^(ud) 01 ud) and 



(3J 



Besides the multiphcative factor of 24 which is discussed in Section formula (|38D establishes 
the relation between the values given in the first line of Table ^ and the values given in Table ^ 
for the following matrix elements: 

{7rrin\ao) = V24:V3{t8\t) = ^ and (ttt/, |ao) = . (39) 

For the coupling of an isodoublet lowest-lying scalar meson to the isodoublet-isoscalar pair, 
we may also select representants. Let us consider the matrix element 



{{us)(f)\us) . 

In this case, the matrix element for the dd contribution vanishes. Consequently, for the isoscalar 
flavour singlets and octets, we end up with 



= sJliABC) and {{us) (l)s\us) = -^^{ABC) , (40) 

which explains the factor 8 in Table | for the quadratic matrix elements under d8 and dl, in the 
sector under isodouhlets. 

From Eq. (^OD, we obtain for the ideally mixed combinations ( pTf ) the relations 



{{us)4)n\us) = y^{BAC) = y^{{us)(f)i\us) and 
= {ABC) = V?((ms)0i , 



which, by the use of the first line of Table ^ and when, moreover, multiplied by the factor v24, 
gives the matrix elements of k ^ Krjn and k, Kr]g, i.e.. 



{Kr]n\K) = V24:^J^{dl\d) = ^\ and {Kr]s\K) = V24^|(rfl |c/) = ys > 

the quadratic sum (= |) of which is found for model B86 in Table ^. 

Now, according to formula (p!?]), for the couplings of lowest-lying scalar isoscalars to meson 
pairs, the normalisation may be not unity. Let us begin with the coupling to a pair of isotriplets, 
e.g. 

((W) {du) 10) . 

When (f) represents a flavour- singlet isoscalar, then, also using formulae ([TTp and (33), one finds 



ud) {du)\(l)i) = ((ud) {du)\^(uu + dd + ss)) = \/||^^^^ + ^^^} = vf^^^C*) 
whereas, when (p represents a flavour-octet isoscalar, it follows that 

[ud) {du) \(j)s) = ^{{ABC) + {BAC)} = ^{ABC) . 



Indeed, in Table | the quadratic matrix elements under tt in the sector for the flavour-octet 
isoscalars are the same as in the sector for the flavour-singlet isoscalars. 



For the ideally mixed combination (|37[), also applying formulae ( ]T7| ) and (0), one obtains 



'{ud){du)\<Pr,) = ^{^^ + ^^} = Vf^^^^) = V1(H Wl</>i> • (41) 



When, moreover, multiplied by the factor v24, formula (41) establishes the relation between the 



first line of Table H and the matrix element of ?7„ tttt, according to 



{7nr\r]n) = v/24J|(tt|0i) - -'^ 



5 ' 



as is found for model B86 in Table ^ 

Next, let us also study the coupling of isoscalars to a pair of isodoublets, e.g. 

{{us) (su) \(f)) . 

When represents a flavour- singlet isoscalar, then, again through the use of formulae ([T7| ) 
and we have 

{{us){su)\^,) = = . 

Similarly, for a flavour-octet isoscalar we find 

{{urs){su)\(l>s) = Jl{{ABC)-2{BAC)} = - Jl{ABC) , 



which result explains the factor ^ in Table ^ between the quadratic matrix elements under dd in 
the sectors for the fiavour-octet isoscalars and for the fiavour-singlet isoscalars. 

For the ideally mixed isoscalars defined in formula (|37|) , using once again formula (^^, we 
obtain 



{{us) {su) = = ^U^BC) = ^l^{{us){su)\<P,) and 



3 

{{u-s) {su) \<f>,) = Vi^^ = ^'ii^BC) = ^l{{u-s){su)\<f>,) 

V 3 

which nicely explains the values for the matrix elements of {rin/^s) ~* KK, i.e.. 



-{BAC) _ ,3 



KK\r]n) = V2Aj^{dd\(f)i) = J\ and {KK\r]s) = V24V|(rfc/ |0i ) 



given for model B86 in Table |^. 

Similar straightforward calculations lead to the matrix elements {rjnTin and {rjsTjs \ris 
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Table 1: Particle identification used in this paper. 
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Table 2: Coupling constants for the decay processes of pseudoscalar mesons into meson pairs. 
The interpretation of the content of the table is explained in the text. 
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Table 3: Coupling constants for the decay processes of vector mesons into meson pairs. The 
interpretation of the content of the table is explained in the text. 
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1 

48 


1010 


1010 


001 


1 

216 


1 
324 


1 

162 


1 

72 


1 
144 


1 

1296 


1 
162 


1 

72 


1 
216 


1 

648 


1 

648 


1 
162 


1 

72 


1 

216 


1 
162 


1 
648 


1 
648 


1 

72 


1010 


1011 


000 


1 

432 


1 
648 


1 

324 


1 

144 


1 
288 


1 

2592 


1 
324 


1 
144 


1 
432 


1 

1296 


1 

1296 


1 
324 


1 
144 


1 

432 


1 

324 


1 

1296 


1 

1296 


1 

144 


1010 


1210 


000 


5 
108 


5 
162 


5 
81 


5 

36 


5 

72 


5 
648 


5 

81 


5 

36 


5 
108 


5 
324 


5 
324 


5 

81 


5 

36 


5 
108 


5 

81 


5 
324 


5 
324 


5 

36 


1100 


1100 


000 


1 

144 


1 
216 


1 

108 


1 

48 


1 

96 


1 

864 


1 

108 


1 

48 


1 
144 


1 
432 


1 

432 


1 
108 


1 

48 


1 
144 


1 
108 


1 
432 


1 
432 


1 

48 


0110 


0110 


000 


1 

48 


1 

72 


1 

36 


1 

16 


1 

32 


1 

288 


1 

36 


1 

16 


1 

48 


1 
144 


1 
144 


1 

36 


1 

16 


1 

48 


1 

36 


1 
144 


1 
144 


1 

16 


1110 


1110 


000 


1 

36 


1 

54 


1 

27 


1 

12 


1 

24 


1 

216 


1 

27 


1 

12 


1 

36 


1 

108 


1 

108 


1 

27 


1 

12 


1 

36 


1 

27 


1 
108 


1 
108 


1 

12 


0000 


1110 


110 


1 

18 


1 

27 


2 

27 


1 
6 


1 

12 


1 

ICS 


2 

27 


1 

6 


1 

18 


1 

54 


1 

54 


2 

27 


1 
6 


1 

18 


2 

27 


1 

54 


1 

54 


1 
6 


1010 


1100 


110 


1 

18 


1 

27 


2 

27 


1 

6 


1 

12 


1 
108 


2 

27 


1 
6 


1 

18 


1 

54 


1 

54 


2 

27 


1 
6 


1 

18 


2 

27 


1 

54 


1 

54 


1 
6 


1010 


1010 


220 


5 

54 


5 

81 


10 
81 


5 

18 


5 

36 


5 
324 


10 
81 


5 

18 


5 

54 


5 
162 


5 
162 


10 
81 


5 

18 


5 

64 


10 
81 


5 
162 


5 
162 


5 

18 



Table 4: Coupling constants for the decay processes of scalar mesons into meson pairs. The 
interpretation of the content of the table is explained in the text. 



initial meson 


decay 

E 


products 

17 


ao or 6 
Kq or K 
fo or e/ S nn 
fo or e/ S ss 


1 (KK) + 1 (7rr7„) 

3 (^^) + 1 (;^;^) + 1 (^^^^^) 
1 (iTiT) + 1 (r^.r^,) 


^ (KK) and ^ (sum vrr^'s) 

2 (-f^Tr) and ^ (sum Krj^s) 

1 (tttt), (KK) and ^ (sum ?7r7's) 

^ (KK) and tj (sum r^r^'s) 



Table 5: Quadratic coupling constants for the decay process of a scalar meson into a pair of 
pseudoscalar mesons. 



